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SIMPLIFIED METHOD FOR CALCULATING THERMAL CONDUCTANCE 
OF ROUGH, NOMINALLY FLAT SURFACES IN HIGH VACUUM 
by Daniel J. McKinzie, Jr. 

Lewis Research Center 

SUMMARY 

A simplified approximate method for calculating the coefficient of thermal contact 
conductance has been derived. The method applies to the case of two rough, nominally 
flat, outgassed surfaces in contact and under load in a high vacuum environment. The 
asperities of the surfaces are assumed to undergo plastic deformation. The method 
does not require the use of a computer program or the calculation of any abstract quan- 
tities. The surface parameters are obtained directly through the use of a surface ana- 
lyzer and a microhardness tester. The derivation is based on the so-called "button" 
model of contact and is done by applying the basic steady- state Fourier heat conduction 
law written at the contact plane. An effective zone of the thermal disturbance is as- 
sumed to exist in the vicinity of the interface. The minimum length of this zone is as- 
sumed to be equal to the summation of the average of the maximum roughness heights of 
the two surfaces in contact. 

The method is applied to an experiment involving two 2. 54-centimeter-diameter 
tungsten cylinders. The cylinders are held in contact and under axial load in a vacuum 
environment of approximately 10 ° torr. The conductance was experimentally obtained 
as a function of the mean interface temperature at various constant contact pressures. 
However, the mathematical form of the experimental data was changed and plotted so 
that the effective area ratio of contact against contact pressure at constant interface 
temperature could be analyzed. This was done to demonstrate the change in the contact 
area with pressure that occurs during the history of such a contact. The agreement be- 
tween the experimental results and theory is considered good. 

A new technique for determining the variation of the specimen's contacting surface 
effective microhardness against temperature is also presented. 



INTRODUCTION 


The prediction of the coefficient of thermal contact conductance across metallic 
joints under load has received much attention during the past two decades, as indicated 
by references 1 to 17. Of particular interest is the work presented by Fenech and 
Roshsenow (ref. 2), Henry and Fenech (ref. 11), and Mikic and Rohsenow (ref. 15). 

Fenech and Rohsenow (ref. 2) have considered the more general problem of deriving 
an expression for the thermal contact conductance of rough, nominally flat surfaces in 
contact as a function of the ambient fluid conductivity. As a limiting case, the authors 
considered the surfaces outgassed to a high vacuum. They chose the thermal conduc- 
tivities of the specimen and the fluid filling the voids, the real area in contact, the num- 
ber of contact points per unit cross-sectional area, and the volume average thickness of 
the void gaps as independent variables. These parameters have the shortcoming that 
they are not accessible to direct measurement but are only calculated or inferred from 
other measurements. Henry and Fenech (ref. 11) present a method for obtaining analog 
voltages of surface profiles of the contacting surfaces. They also present a general pur- 
pose analog computer program to determine the values of the geometric parameters of 
the contact interface (proposed by Fenech and Rohsenow) as a function of contact pres- 
sure. Thus, Henry and Fenech combine their computer program with the analysis of 
Fenech and Rohsenow. Mikic and Rohsenow (ref. 15) have treated the surfaces of con- 
tact as random processes with Gaussian distributions of height. They have considered 
four cases in particular. The solution they present for the case of rough, nominally flat 
surfaces in contact in a high vacuum environment is in a more simple form than that 
presented by Fenech and Rohsenow. The method of Mikic and Rohsenow requires the 
limited use of a computer in addition to the surface analyzer and surface microhardness 
tester. Though these analyses appear to have been experimentally confirmed by their 
authors, other investigators have found that the computerized programs outlined in ref- 
erences 2, 11, and 15 have required greater effort and cost than simply obtaining test 
data directly from sample specimens. 

The investigators of references 2, 3, 5 to 7, 10 to 15, and 17 claim good experi- 
mental and theoretical agreement from their own experiment. The experimenters of 
references 1, 8, 10, and 13 have had difficulties in applying the theories of others to 
their test data. In part, this lack of agreement is believed due to the following: 

(1) There is insufficient knowledge of the microscopic and macroscopic structure of 
the specimen's contacting surfaces, where this knowledge is required to evaluate com- 
plex parameters appearing in the analyses. 

(2) There exists a lack of similarity between theoretical and experimental models. 

(3) Little or no consideration is given to the discrimination between elastic or 
plastic deformation phenomena that take place between the surface asperities of the test 
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specimen. If the conditions that bring about elastic or plastic deformation are not under- 
stood, their occurrence may become confused. This can happen, for example, if the ex- 
perimenter inadvertently loads the asperities of his specimen well into their plastic de- 
formation regime and then unloads them before gathering test data. The initial elastic 
phenomenon is obscured in this case because of the overshoot into the plastic regime 
where the microstructure of the surface asperities and properties of the surface mate- 
rial have changed irreversibly from their unloaded condition. Such an event may take 
place in the experimenter's attempt to aline the test specimen interfaces. Thus, the 
initial elastic and plastic deformation phenomena would not be observed. An additional 
elasto-plastic problem may occur if the experimenter allows a change to take place in 
the mean interface temperature of his specimen during a test without considering the 
effect this change has on the hardness properties of the material. 

Because of these problems and the complications introduced by the need for the 
computerized theoretical approaches of references 2, 11, and 15, the purposes of the 
present work are as follows : 

(1) Reemphasis is made of the description of the elastic and plastic deformation of 
surface asperities on nonmoving surfaces in contact and under load as presented by 
Bowden and Tabor (no thermal effects included) in reference 18. This reemphasis is 
done to clarify the proper sequence of loading the interface so that meaningful experi- 
mental data can be obtained. 

(2) Reemphasis is made of the importance of knowing the history of the contact joint. 

(3) A new approximate method is presented for predicting the coefficient of thermal 
contact conductance for rough, nominally flat, outgassed surfaces undergoing plastic de- 
formation in a high vacuum environment. This new method, which does not require a 
computer evaluation of the coefficient of thermal contact conductance, was derived by 
applying the basic steady- state Fourier heat conduction law to the so-called "button” 
model of contact. 

(4) The new approximate method is applied to the interface of a tungsten-tungsten 
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specimen tested in a high vacuum environment of approximately 10 torr, with ground 
surfaces of contact and approximately parallel lays. 

(5) A new method is presented for determining the variation of the effective micro- 
hardness with temperature of the interface surface asperities. 


MECHANICAL CONTACT PHENOMENON 

The following discussion of the mechanical contact phenomenon, which is presented 
in greater detail in reference 18, consists of a description of the mechanical aspects of 
the loading, work hardening, and unloading of a surface asperity. This discussion is 
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included because it describes the reversible and irreversible events occurring during 
the making of a contact. It also demonstrates the importance of knowing the history of 
the contact event. 


Loading of an Asperity 


Reference 18 considers an idealized case where the tip of an asperity is assumed to 
be smooth and spherical in shape. A single asperity is assumed to rest on a softer 
metal, and, in the region of contact, the surface of the softer metal is considered to be 
a plane (see fig. 1(a)). Consider what happens as the load on the asperity is increased 
in three steps. 

First, if the surface asperity is pressed into the plane with a load W', it at first 
deforms elastically according to Hertz's classical equations. (All symbols are defined 
in appendix A. ) At this stage, the planar projection of the area of contact A’ (a' = v a 2 ) 
is proportional to W' / , since from Hertz's equation written for a flat surface being 
penetrated by a single spherical asperity 


a - 1.1 



1/3 


( 1 ) 


where a is the radius of the contact indent (see fig. 1(b)), r is the initial radius of 
curvature of the asperity, and and Eg are Young's moduli for the two surfaces, 
respectively. In this range the deformations are reversible. 

Second, as the load is increased, the mean pressure p m (i. e. , p m = W' /A' ) over 
the contact area of the asperity increases until it reaches such a value that at a critical 
point z (see fig. 1(c)) within the softer metal the elastic limit is exceeded. From refer- 
ence 18, the Hertzian analysis indicates this region is situated at a point about 0. 6a be- 
low the center of the area of contact. At this stage the metal reaches a state of plasticity 
and yields irreversibly. The metal outside the vicinity of point z has not yet reached 
the conditions for plasticity, so its deformation is still essentially elastic. 

Third, as the load is increased further, the area of contact A^ and the mean pres- 
sure p m rise in a manner which deviates increasingly from the elastic relation. The 
region of plasticity at point z below the center of the area of contact grows rapidly, and 
a stage is soon reached at which the whole of the metal around the region of contact is 
flowing plastically (see fig. 1(d)). At this condition, the theoretical works of Hencky and 
of Ishlinsky mentioned in reference 18 show that 

P m = cY (2) 
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n 


where Y is the elastic limit of the softer metal as found in pure tension, and c has 
a value of approximately 3. 


Work Hardening of an Asperity 

Bowden and Tabor (ref. 18) consider the effects of work hardening in an attempt to 
describe the change of A' c with increased p m> They point out that all metals experience 
work hardening during increased loading. As before, they assume a hard smooth spher- 
ical asperity penetrates the surface of a smooth softer metal. Their assumed sequence 
of events is now described. 

As a load is applied, the formation of an indentation produces an increase in the 
elastic limit Y of the surface metal. It follows from this result that when metal sur- 
faces rest on one another the peaks of their asperities readily deform plastically with 
the mean pressure given by equation (2). The parameter Y in equation (2) is some 
"representative" measure of the elastic limit of the deformed metal at the tip of the as- 
perities, and c depends on the shape and size of the surface irregularities. The coeffi- 
cient c has a value of about 3 for hemispherical asperities and for conical and pyra- 
midal asperities over a wide range of angles. 

With conical and pyramidal asperities the yield pressure is independent of the 
amount of deformation that has occurred. Thus, the area of real contact A^ is directly 
proportional to the load W'. With hemispherical asperities this is only true for metals 
whose surfaces are highly work hardened. 

For all the surfaces considered herein the surface asperities are assumed to be 
highly work hardened by the grinding process used to prepare them. Therefore, for all 
types and shapes of surface irregularities, the real area of contact A^, is considered 
proportional to the load in the plastic regime. Thus, 

A’ ocW’ (3a) 

or 

A* 

— «: p’ (3b) 

A 

where 



A 


and A is the total cross-sectional area of the surface considered. 
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Unloading of an Asperity 


Next, the unloading of the surface asperity is considered. From reference 18, in 
the range of initial elastic deformation, the process of deformation is completely re- 
versible. Hence, a removal of the load enables the contacting surfaces to return to 
their original configurations. 

For the case of initial plastic deformation of a single asperity, consider a hard 
sphere of radius a penetrating a softer metal under a load W’ Q to form an indentation 
of diameter d’ (circular planar projection of the indent). Suppose the load on the 
sphere is removed and then replaced. The indentation should now be in exactly the same 
condition as when it was originally plastically formed. Thus, any change which has oc- 
curred in the configuration of the surface between removing and replacing the indenter 
over the same loading range must be reversible; that is, it must be elastic. Conse- 
quently, the recovered (unloaded) indentation has a radius of curvature larger than that 
of the spherical indenter. At any stage in this process, if the load is reduced, there is 
a relaxation of elastic stresses, and the area over which the surfaces touch is given by 
the laws of elastic deformation. It should be pointed out that this statement does not 
imply that equation (1) still holds. Specifically, equation (1) is written for a single as- 
perity penetrating a flat plate and not for a single asperity penetrating a surface which 
is concave, as would be the case here. It does, however, imply that an indentation di- 
ameter d' is reached at the stage where the load returns to the value W' Q and the 
elastic deformation has reached its end. The whole of the surface is in a stage of in- 
cipient plasticity, and, if the load is increased by a small amount, further plastic de- 
formation occurs and the indentation increases in size. 

The process of loading and unloading an asperity is described graphically in fig- 
ure 2, which is similar to a figure given in reference 18. Figure 2 indicates that, as 
the load on a single asperity is increased, an initial elastic deformation of the area of 
contact occurs; this is predicted by equation (1). At some larger load the onset of plas- 
tic deformation takes place, and the variation of the area of contact A^ with load is 
then predicted by the approximate equation (eq. (2)) for plastic deformation. If the as- 
perity is then unloaded at some higher load W' Q and area of contact A' c , the variation 
of the area of contact against load is indicated by the dashed curve labeled reversible 
elastic deformation. If, further, the load on the asperity is removed, an offset in the 
area of contact A' « due to plastic deformation will occur, as shown in figure 2. 

C, 1 


THEORETICAL ANALYSIS 


The object of this analysis is to determine an approximate parametric relation for 
the coefficient of thermal contact conductance written for two rough, nominally flat, 
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outgassed surfaces in contact in a high vacuum environment, where the asperities on 
these surfaces are undergoing plastic deformation. The analysis is based on the steady- 
state Fourier one- dimensional heat conduction law. It is applied to the so-called 
"button" model of contact described in reference 2. A discussion of the experimental 
determination of the quantities appearing in the derivation is included. Also included is 
a parametric relation between the coefficient of the thermal contact conductance and the 
contact pressure for the unloading of an interface. 


Coefficient of Thermal Contact Conductance for Plastic Regime 

Figure 3 is a view of the assumed contact interface showing a contact region labeled 
the effective zone of the thermal disturbance. The length of this zone is assumed to in- 
clude the sum of the largest asperity roughness heights of both portions of the specimen. 
This height is used as the best measurable length approximating the physical length of 
the zone of thermal disturbance. Because this length may be smaller than the actual 
length over which the thermal discontinuity at the interface takes place, it is referred to 
here as the effective zone of the thermal disturbance. 

Assumptions . - The following assumptions must be applicable to the test specimen 
before the parametric relations developed in this work may be used: 

(1) The ends of the cylindrical specimen are nominally flat (as defined in ref. 6). 

The surface features are measured in microinches and defined in reference 19. 

(2) Plastic deformation of the asperities in contact is taking place. 

(3) The surfaces are clean and outgassed to a high vacuum environment. 

(4) The contacting cylinders are made of the same material. 

(5) Effects of creep are negligible across the interface of contact. 

(6) Effects of radiation are negligible across the interface of contact. 

(7) The outside walls of the specimen are adiabatic. 

(8) The axial length of the effective zone of the thermal disturbance 6 is very small 
compared to the specimen’s diameter. 

The following assumptions apply to the theoretical model which is used to approxi- 
mate the specimen. 

(1) The specimen's contacting surface asperities are each replaced by the so-called 
"button" model of contact. 

(2) Each button contact is thermally insulated from its neighbor. 

(3) The number of button contacts is uniformly distributed over the interface of con- 
tact. 

(4) The axial temperature distribution in each cylinder of a button contact is the 
same as every other one. Thus, their circular areas of contact may be summed. The 
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resulting single area of contact is shown in figure 4 (modified from ref. .2) and shall 
hereafter be referred to as "the button contact. " 

(5) The face of the button contact is in a state of plastic deformation with increasing 
load. 

(6) No thermal gradient exists in the radial direction of the cylindrical button con- 
tact. 

(7) The axial temperature gradient is constant and remains steady in time throughout 
the cylindrical volume of the button contact. 

(8) The length of the effective zone of the thermal disturbance is very small and is 
equal to the sum of the lengths of the button contacts of the upper and lower portions of 
the model. It is therefore assumed that only within the effective zone of the thermal dis- 
turbance does the thermal discontinuity at the interface exist and that the thermal gradi- 
ent in this zone is constant and steady in time. This concept is unique to the present 
analysis. 

(9) The interface of contact between the button contacts is located at a point one-half 
the axial length of 5 (see fig. 5). This assumption is made though the lengths of the 
button contacts 6^ and shown in figure 4 are not necessarily equal. This relocation 
does not introduce an error into the theoretical analysis because the length of <5 is as- 
sumed to be very small. Thus the change in the mean temperature of the interface of 
contact (fig. 4) is negligible. 

(10) Outside the effective zone of the thermal disturbance the temperature distribution 
is steady in time throughout the cylindrical model. 

Analysis . - When these assumptions are applied to the model, the axial temperature 
distribution in the vicinity of the interface is as shown in figure 5. The basic steady- 
state one- dimensional Fourier heat conduction law written for the effective zone of the 
thermal disturbance 6 (figs. 4 and 5) is given by 

-9_=-k^ (4) 

A c 

where A is the sum of the individual cross-sectional areas of the idealized individual 
c 

"button" contacts, and k is the coefficient of thermal conductivity. This coefficient is 

evaluated at a temperature equal to (t. . + t, 0 )/2 (see fig. 5), and dt/dx is the axial 

i, 1 1 > A 

temperature gradient across the disturbance zone 5. 

If equation (4) is applied to the effective zone of the thermal disturbance (see fig. 5), 
the heat-transfer rate per unit area at the interface of contact may be approximated by 

-k At 

Q = £ (5) 

A 5 
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where from figure 5 at x = -x^ 


and at x = Xg 


t = t (D 
t = t d) 


Thus, 


At c = *(!) " *0) 


( 6 ) 

(V 

( 8 ) 


and equation (5) becomes 


Q . ~ t (£)) 


A 


c 


6 


(9) 


The coefficient of thermal contact conductance h is defined in the equation 


— = -h At (10) 

A 


where A, the cross-sectional area of the specimen, is expressed by 

A = ttR 2 (11) 

If the temperature gradients in the specimen outside the effective zone of the thermal 
disturbance are linearly extrapolated to the interface x = 0, as shown in figure 5, then 
At at the interface is given by 


At = ti,2 -ti,i 


( 12 ) 


Substituting equation (12) into equation (10) gives 


Q 


--Wi.a-*!,!) 


(13) 


Equating the heat-transfer rate Q across the interface, equations (9) and (13), yields 


9 


(14) 


kA. 


ft® ' t ©) _ 


= Ah(t. 9 -t, J 


6 i,2 1,1' 

From figures 4 and 5 and equations (8) and (12), 

l^t c | S: | At | 

Since 6 is very small, it is reasonable to assume 


l*t c ] = I At | 


Substituting equations (8) and (12) into equation (16) yields 


‘ t <D ~ \ 2 " \ 1 


(15) 


(16) 


(17) 


Equation (14) thus reduces to 



k A 


From plastic theory for nominally flat surfaces in contact (ref. 2), 


A I 


(18) 


(19) 


where P is the pressure resulting from the applied load W acting over the total nomi- 
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nally flat surface tfR , and I is the average of the random microhardness values of the 
softer of the two interface surfaces. Substituting equation (19) into equation (18) results 
in the desired relation 


h - = — (20) 

k I 

Unloading the Contacting Surfaces 

Consider the case of unloading the contacting surfaces from some plastically de- 
forming pressure P Q at a constant mean interface temperature where P Q = (w q j ?rR^) . 
Equation (1) can be applied to the summation of the individual contact areas. Therefore, 
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( 21 ) 


A ocw 2 / 3 
c 

Multiplying equation (21) by 1/A yields 

A c „ W 2 / 3 

A A 


( 22 ) 


o 

where A is rrR . If some upper load W Q is set as a limit, the area of contact A c 
has the same magnitude for the elastic and plastic regime. Substituting equation (18) 
into proportionality (22) at load W Q results in 


h 6 W 2 / 3 

_2_ cc 9 

k A 


(23) 


or 


V 

k 


p 2/3 

o 

a 1/3 


(24) 


In the elastic regime at a constant mean interface temperature assume 


£ = C (25) 

k 


where C is a constant. If now the interface is unloaded from P„ to P, 

o 1 



(26) 


where A Q = A^ 


and C Q = Cj. 


Equation (26) thus reduces to 



(27) 
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Substituting equation (20) into equation (27) yields 


h. = p2/ 3 pV 3 JL 
1 1 ° 61 


( 28 ) 


or substituting equation (18) into equation (27) yields 

I A _ 



(29) 


Thus equation (28) provides a method for calculating the coefficient of thermal con- 
tact conductance for any point below the upper limit of P Q and h Q . This is true if the 
mean interface temperature of the contacting surfaces and the heights of the plastically 
deformed asperities remain constant or very nearly constant during the unloading pro- 
cesses. 


EXPERIMENTAL APPARATUS AND PROCEDURE 
Test Apparatus 

Sketches of the bell jar, its contents, and its location in the test setup are shown in 
figures 6 and 7. Simply, the test is performed in a bell jar that contains an induction 
heating coil and a column composed of the test specimen, a heat meter, and heat sink as 
shown in figure 6. A pneumatic loading cylinder is located on top of the bell jar as 
shown in figure 7. The loading cylinder is used to vary the pressure at the contact inter- 
face of the test specimen. Electrical currents are induced in the extreme upper portion 
of the test specimen. The heat generated by these currents is transferred across the 
test interface, through the heat meter, and finally to the water-cooled heat sink. The 
test specimen and heat meter are instrumented with thermocouples. The millivolt level 
output of the thermocouples is read out on a potentiometer. The bell jar operates at a 
vacuum condition in the 10 -D torr range. A detailed discussion of the test apparatus and 
specimen is given in appendix B. 


Experimental Determination of h 

The experimental method used to determine the coefficient of thermal contact con- 
ductance h of a tungsten-tungsten interface is given in this section. 
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Test procedure . - The coefficient of thermal contact conductance of the tungsten- 
tungsten specimen was experimentally determined using equation (13) and the test appa- 
ratus described in appendix B. Heat was transferred from a source across the tungsten- 
tungsten interface of contact and then to a heat sink (see fig. 6). The rate of heat trans- 
fer across the interface was determined with the aid of an oxygen-free high- conductivity 
(OFHC) copper heat meter. The test specimen was designed so that a one-dimensional 
linear heat-transfer analysis would apply. Thus the hypothetical At of equation (10) 
was experimentally determined by a linear extrapolation of the thermal gradients from 
the upper and lower portion of the specimen to the interface (see fig. 5). 

Before proceeding it should be stated that the sequence of changes made in pressure 
and temperature were not ideally made throughout the test. For the first 22 runs, out of 
a total of 72, the data were obtained by correctly varying the test conditions. As noted 
in figure 8, these data were obtained by correctly increasing the pressure between suc- 
cessive testing periods in which the interface temperature was changed in a monotonic 
increasing manner. Thus changes in the area of contact due to elastic or plastic effects 
taking place with a rise in interface pressure and/or temperature were distinctly ob- 
tained. The data of figure 8 were subsequently cross plotted in figure 9 with the aid of 
equation (18). Figures 10 and 11 present the data of figure 9 in varying degrees of re- 
finement. 

Three power failures occurred during the test resulting in four subgroups of data, 
as noted in figure 11. However, the contact pressure remained constant during each of 
them. The first two interruptions occurred after runs 12 and 32 when the interface tem- 
peratures were approximately 331 and 382 K, respectively. And the third occurred after 
run 36 when the interface temperature was 529 K. After each interruption the test spec- 
imen cooled to room temperature before power was reapplied. Thus differences in lat- 
eral contractions resulted in sliding among the asperities. 

The shortest period of time needed to establish steady- state conditions after a 

change was made in pressure was approximately lw hours. The average time required 

1 ^ 

was 2~2 hours. Changes in interface temperature were made in a minimum time of 

3 hours, with the maximum being 48 hours. These intervals of time were considered 
adequate based on previous experience. However, better steady-state data would have 
been obtained with finer control of the power. 

The interface surfaces were cleaned with acetone before installation in the test fa- 
cility. In future tests the necessity for using some standard cleaning procedure seems 
apparent, because of the possibility of film deposition from such a solvent as acetone. 

Data reduction . - The upper piece of the test specimen (fig. 6) was instrumented 
with six special type Chromel-Alumel 30 gage wire thermocouples. The lower piece 
was instrumented with five thermocouples. Six more were located in the OFHC copper 
heat meter at three axial stations. The value of the thermal conductivity used to com- 
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pute the heat transfer across the heat meter was obtained from data presented in refer- 
ence 20 and is shown in figure 12. Figure 13 demonstrates a typical best line fit drawn 
through the experimental thermocouple data. The linear extrapolation of these data was 
used to determine the interface temperature difference. Two typical runs, runs 13 
and 23, are shown. The thermocouple readings at the 3.2 millimeter locations on either 
side of the specimen’s interface were omitted from consideration. This was done be- 
cause of the constrictive resistance effects (ref. 15) taking place in the vicinity of the 
interface. A further simplification was made in the data reduction by determining the 
thermal gradient in each portion of the specimen from two thermocouple readings in the 
same manner as that used by Fried (ref. 2). The gradients determined by this technique 
were continually checked by plots such as that of figure 13. This simplified procedure 
used the same thermocouples throughout the tests. Thus, an added degree of consistent 
variation in the data resulted with varying input heating rates. The error in the temper- 
ature difference across the interface introduced by this method was generally less than 
1. 1 K from a best line fit through all the thermocouple data. However, the manufactur- 
er's accuracy tolerance of ±1. 1 K must apply in an error analysis because a calibration 
of the thermocouples on the specimen was not made. It should be mentioned that the 
large temperature difference, At, obtained across the interface during the majority of 
test runs helped to minimize the error introduced by using Fried's first-order fit. The 
smallest value of At was 9 K which occurred at a mean interface temperature of 331 K 
and a pressure of approximately 5. 31x10° newtons per square meter. This occurred for 
1 run out of the total of 72. The next smallest At was 18. 3 K which also occurred only 

once. The largest value of At, 250 K, occurred at a mean interface temperature of 
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666 K and a pressure of approximately 4. 83x10 newtons per square meter. 

Theoretical Determination of h 

The value of h was determined from values of 6, k, and I together with equa- 
tion (20). The quantities 6 and I were evaluated from the test specimen. The varia- 
tion with temperature of tungsten's thermal conductivity k, which was estimated from 
the literature (ref. 21), is presented in figure 14. 

Measurement of the effective zone of the thermal disturba nce 5. - A Brush Surfan- 
analyzer was used to determine the magnitude of 6. It consists of a control unit, a 
stylus 1200 drive unit, a 50-milligram stylus probe, and a recorder unit. This appara- 
tus records the surface roughness (width and height), profile, waviness, and arithmetic 
average height. All of these quantities are defined in reference 19. 

During the approach of two rough, nominally flat surfaces the maximum roughness 
height would make contact first. Therefore, it is appropriate to determine the average 
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value of the maximum roughness heights (upper) and (lower) of the test specimen 
before the test is made. Since it is not possible to obtain surface analyzer traces which 
include all the points of contact on a surface, a representative number of traces must be 
taken. From each of these traces, a single value of the maximum roughness height is 
obtained. These values are then averaged to obtain the value of the roughness height 
used in the calculations. The value of 6, where 5 = 5^ + Sg, shall hereafter be re- 
ferred to as the effective value of the thermal disturbance height. For the test reported 
here the magnitudes of and 5g (5.35 and 5.51 jum), respectively, were determined 
by averaging the maximum roughness heights obtained from six traversals at different 
positions on each face of the specimen perpendicular to their ground lay. 

In the test reported here the magnitudes of and 6g were determined after the 
test data were obtained. This occurred because the instrumentation required to deter- 
mine and Sg was not available before the test was performed. Ideally, the maxi- 
mum roughness height should be determined from the specimen before the conductance 
test is started. If, however, the test is completed before the traces are obtained, ref- 
erences 2 and 22 indicate that the roughness heights will be only slightly smaller. This 
is only true, however, if full plastic flow of the cylindrical specimen is far from being 
approached. The loads applied and of interest in the test considered here, and most 
other contact conductance investigations, satisfy this condition. In references 2 and 22 
the specimen’s roughness height did not decrease further after its initial plastic defor- 
mation even though much larger loads were applied. 

As shown in figures 8 and 11, the arithmetic average roughness heights for the 
upper and lower surfaces were 0.538 and 0.943 micrometer, respectively. 


Measurement of the Microhardness 

As pointed out in reference 3, the microhardness of cold- rolled nickel and mild 
steel contact surfaces was up to five times larger than their Meyers hardness value de- 
termined from a 1.905-centimeter-diameter cylinder. Specifically, the microhardness 

4 

increased as the load on the microhardness indenter decreased from 4.44x10 to 0.444 
newton. Reference 3 indicates that Laming (ref. 23) observed an apparent increase of 
the microhardness with decrease in the indenter loading for steel-brass and steel- 
aluminum test specimens. Reference 11 also indicates an increase in the microhardness 
of a surface as the indenter load is decreased. 

The test reported here has also indicated an increase in the value of the contact sur- 
face microhardness as the indenter load is decreased. The increase in microhardness 
was obtained with a Knoop indenter at room temperature and is shown in figure 15. The 
choice of the Knoop indenter used in this test was based on the advantages of its use that 
were mentioned in reference 24. The increase in surface microhardness is attributed to 
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surface work hardening caused by the method of surface preparation (in this test, ground). 
This increase has proved to be of great importance. As pointed out in reference 11, the 
correct value of the microhardness is that value obtained when the hardness indenter 
penetrates to the same depth that an asperity penetrates the specimen's surface. How- 
ever, the depth of asperity penetration is a difficult measurement to make. Because the 
penetration depth of the Knoop indenter is geometrically related to the length of its major 
and minor diagonals, these observed indents may also be used as a measure of the depth 
of indenter penetration. This is done by comparing at the same magnification a metal- 
lograph photograph of the micro hardness tester indents at various loads with one of the 
asperity indents. The larger diagonal of the Knoop indenter is best used as an indicator 
of the penetration depth. This is because the included angle of the indenter which forms 

1 ° 

the major diagonal is approximately 172^ , which is more closely equal to the included 
angle of most characteristic asperities (160° to 164°, ref. 6). This range of included 
angles has also been confirmed by our experience. 

Although this technique of making the indenter penetration comparable to the asperity 
is desirable, it could not be used entirely in the test reported here. This was because 
the Tukon microhardness tester used had a minimum Knoop indenter load of 0. 245 new- 
ton. At this load the larger diagonal indicated the indenter penetration depth varied from 
about two to five times greater than the average asperity indent. However, the value of 
the hardness determined with this load was used in the analysis. The variation of the 
microhardness with indenter load is presented in figure 15. 

Effect of Temperature on Hardness 

The tensile strength, thus the hardness, and recrystallization of the tungsten speci- 
men are functions of temperature (ref. 25). Therefore, it is important to include the 
effect of temperature on the hardness in the calculation of the theoretical value of the 
coefficient of thermal contact conductance. The variation of the specimen's surface 
microhardness with temperature was not obtained from microhardness data. However, 
it is suggested that this variation may be determined by applying the subsequent proce- 
dure. If it is assumed that (1) the specimen's interface asperities are kept under a con- 
stant pressure as the mean interface temperature is raised and (2) the average value of 
the maximum roughness height of each portion of the specimen remains constant during 
the change in temperature, then equation (20) reduces to 
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(30) 



where the subscripts 1 and 2 refer to different temperatures at which data is taken. 
Thus, the variation of k/I with temperature can be obtained as a function of the con- 
ductance. If the thermal conductivity of the specimen is known as a function of temper- 
ature, the variation of an effective microhardness with temperature can be obtained. 
Thus, 


(^\ = — — (31) 

VVeff k 2 h l 

The procedure is best done after the variation of the coefficient of thermal contact con- 
ductance with pressure, at constant interface temperature, has been obtained. And the 
data are best obtained at test conditions slightly higher in pressure and temperature than 
the largest values of these variables used to obtain the conductance data. This is be- 
cause equation (31) is written for the plastic deformation regime. Also it is important 
to keep the specimen's interface temperature well below the temperature of recrystal- 
lization if the specimen is not 100 percent recrystallized. This restriction is important 
because the change in grain boundaries that occurs during recrystallization may have 
adverse effects on the transmission of heat across them. 

The technique was applied to the data of runs 7 to 11 indicated in figure 8 at a pres- 
sure of 5. 31x10 newtons per square meter. Substituting the variation of the thermal 
contact conductance and thermal conductivity with temperature into equation (31) results 
in the variation of the specimen's effective surface microhardness with temperature as 

shown in figure 16, curve A. This curve is referenced to the microhardness obtained at 

9 2 

room temperature and shown in figure 15 (8.72x10 N/m ). The values of the effective 
microhardness used to determine the curves and points labeled as plastic theory in fig- 
ure 11 were taken from figure 16, curve A. The data of curve B, taken from refer- 
ence 26, approximate those of curve A as noted from figure 16. 


DISCUSSION OF RESULTS 

The original purpose of the experiment was to show how the coefficient of thermal 
contact conductance for two rough, nominally flat, outgassed surfaces in contact and 
under load varied with changes in the mean interface temperature of the specimen. The 
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contact conductance data were obtained at interface pressures of 4. 83x10®, 2. 55x10®, 
and 5. 31x10 newtons per square meter (fig. 8). At the time the experiment was per- 
formed it was not realized that a monotonic increase of interface pressure and/or tem- 
perature was the ideal procedure for obtaining data for a specimen. This scheme of 
testing is necessary because of the unique elastic and plastic deformation phenomena 
that occur only during the first loading cycle. These phenomena are thus associated 
with the history of changes in loading and/or temperature of the asperities in contact. 
Not all of the test data presented here were obtained by the ideal variation of pressure 
and temperature. However, the data resulting from the nonideal experimental proce- 
dures used clearly show effects which are explained by theoretical considerations. In 
turn, these results indicate the importance of knowing the history of the changes through 
which the asperities go. 

The data, which are shown in figure 8, are cross -plotted in figure 9 with the aid of 
equation (18). Figures 10 and 11 present the data of figure 9 in varying degrees of re- 
finement. The data of figure 8 were replotted because their form was not easily related 
to the description of the contact phenomena presented here and shown graphically in fig- 
ure 2. Therefore, figures 9 to 11 are plots of the effective area ratio of contact against 
contact pressure. These figures indicate that, at a constant interface temperature, de- 
formation ranging from plastic to elastic took place during the loading or unloading, re- 
spectively, of the tungsten-tungsten specimen. This deformation range supports the 
description given in reference 18. 

The data of figure 8 were obtained in a timewise sequence of five distinct groups. 
Each group consists of data obtained by varying the specimen’s interface temperature 
from a lower to an upper limit. The strict maintenance of constant interface pressure 
for each data group was not considered important at the time the data were obtained. 
Thus, some of the data include small pressure and thermal variations. This testing 
technique induced unwanted effects in the data and would not be used at present. It was 
during the first 22 runs (fig. 8) that both the highest pressure and temperature levels 
were imposed on the specimen. Thus, the initial elastic and plastic deformation of the 
test specimen occurred during these runs. 

The first timewise group of data, A to B, shown as the lowest curve in figure 9 has 
a distinctly different slope of the effective area ratio of contact against increasing pres- 
sure than the remaining groups. Also, the magnitude of the slope for the first group 
agrees with plastic theory, while the slope for the remaining groups appear to agree 
with the elastic theory, though this is not clearly discernible in figure 9. 

Figure 10, which is an attempt to present more clearly the conductance data than is 
done in figure 9, is a plot of selected data from figure 8 at the specific interface temper- 
atures of approximately 378, 392, 529, 570, 628, and 666 K. These temperatures were 
chosen because of the availability of data at two or three pressures. The data of fig- 
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ure 10 indicate the slope of the three point pressure data (approximately 4. 83x10^, 

2. 55x10 , and 5. 31x10° N/m ) is generally not constant with pressure. The slope of the 
data in the supposedly elastic region appears to vary from the anticipated theoretical 
value of 2/3, while the slope of the data between runs 1 and 6 (at an interface tempera- 
ture of approximately 390 K) agrees closely with the value of 1 from plastic theory. 

Thus a somewhat confused picture of just what is happening at the specimen's interface 
of contact with changes in pressure and temperature would result from only considering 
the data presented in figure 10. 

Figure 11, which is a refinement of the data presented in figure 10, presents only 
data obtained after the specimen is believed to have reached thermal equilibrium; in 
addition, possible effects due to several power failures that occurred during the test pe- 
riod are eliminated. Three power failures occurred which resulted in four unique sub- 
groups of data. The power failures occurred at interface temperatures of approximately 
331, 382, and 529 K after runs 12, 32, and 36, respectively. After each, the specimen 
cooled to room temperature while the contact pressure remained constant. The cause of 
these power failures is not known. Because of the likelihood that microscopic lateral 
asperity misalinements took place, due to large interface differential temperature ef- 
fects, it is assumed that the data for each of these groups is independent of the others. 
The curves of figure 11 result by joining like temperature data in each of the subgroups. 
Separating the data in this manner results in slopes for the curves that agree well with 
plastic or elastic theory. 

The variation of the surface microhardness with temperature is necessary to calcu- 
late the effective area ratio of contact for the plastic regime shown in figure 11. This 
variation in microhardness is normally obtained by indenter techniques; however, it 
may be obtained using a new technique which will be described. The specimen’s micro- 
hardness may be calculated as a function of temperature by substituting the experimental 
data of figure 8, the thermal conductivity from figure 14, and a reference value of the 
specimen's microhardness into equation (31). The reference value of the microhardness 
(8. 74x10 N/m ) was obtained by indenter techniques at approximately 294 K and is pre- 
sented in figure 15. The resulting variation of the specimen's effective microhardness 
with temperature (curve A) is presented in figure 16. Also included in figure 16 are 
tungsten hardness data obtained from thin sheet (0.025 cm) highly cold worked high purity 
tungsten having no stress relief (ref. 25) and 0.071-centimeter-diameter drawn high pu- 
rity tungsten wire (ref. 26). These curves were obtained by converting the yield strength 
data presented in these references into hardness data. The conversion was done by 
multiplying the yield strength by the usual coefficient of 3. It should be pointed out that 
the temperature range of the data presented in figure 16 is well below the recrystalliza- 
tion temperature of tungsten and also below the temperature range of stress relief. 
Therefore, the difference in hardness levels is due to the degree of cold working of the 
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tungsten material. It is noted that the slopes of all three curves are approximately 
equal. Thus it appears that curve A represents the hardness variation of the surface 
asperities of the tungsten specimen with temperature. Because of the good agreement 
in magnitude arid slope between the hardness data of curves A and B (ref. 26), the curves 
and points of figure 11 labeled plastic theory were determined by using equation (19) 
where the value of the microhardness was taken from curve A (fig. 16). The agreement 
is considered good between the theoretically calculated and experimentally determined 
values of the effective area ratio of contact against pressure presented in figure 11. It 
is believed that the use of this technique for predicting the variation of the effective 
microhardness of a contact interface is unique to this report. 

As a further check on the plastic theory represented by equation (20), a calculation 
of the thermal contact conductance was made at a pressure of 5. 31x10 newtons per 
square meter and an interface temperature of 294 K (room temperature). The hardness 
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value of 8. 74x10 newtons per square meter was obtained from microhardness indenter 
data (fig. 15), and the conductivity was obtained from figure 14. The value of the calcu- 
lated conductance (9. 11 kW/(m )(K)) is approximately 15 percent larger than that ob- 
tained from the conductance experimental data shown in figure 8 (hj = 7. 93 kW/(m )(K)). 
This experimental value was obtained by extrapolation of the data for runs 6 to 11 to an 
interface temperature of 294 K. These data were used because they are believed to be 
the result of pure plastic deformation of the specimen's interface surfaces of contact. 
Most of the remaining data in figure 8 are believed the result of elastic deformation; 
therefore, equation (20) cannot be used with these data to predict the coefficient of ther- 
mal contact conductance. 

The elastic data shown in figure 11 were obtained at interface pressures below the 
maximum impressed during the initial plastic deformation of the asperities. The mag- 
nitude of the slopes agrees well with the elastic theoretical value of 2/3, as noted in fig- 
ure 11. 

A small error study was performed to indicate the reliability of the theoretical and 
experimental results. The error introduced into the effective area ratio of contact 
varied from approximately -6 percent (below) to 11 percent (above) the theoretical val- 
ues. This is due to the known scatter in the microhardness of the softer of the speci- 
men's surfaces of contact (see fig. 15). The effect of this scatter on the plastic theo- 
retical curve is shown in figure 11 at a mean interface temperature of 390 K. 

No accurate estimate of the error in the effective area ratio of contact for the ex- 
perimental data could be made. This was due to the lack of thermocouple and thermal 
conductivity calibration data of the tungsten-tungsten specimen. The calculation of the 
experimental error was made with the aid of references 20, 21, and 27. When conser- 
vative estimates of the errors are made, the largest error in the effective area ratio of 
contact varied from -18 percent (below) to 10 percent (above) the experimental values. 
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This error, however, applied to only three runs out of the total of 72. The next closest 
estimate of the error resulted in a variation from -11. 3 percent (below) to -0.6 percent 
(below) the experimental values for a total of 10 runs. And the majority of runs had a 
variation from -8.2 percent (below) to -1.6 percent (below) the experimental values. 
Thus, the majority of the data had a mean error of approximately -5 percent (below) the 
experimental values. 

Radiation losses affecting the accurate determination of heat transfer across the 
interface of contact were found to be negligible at all test conditions. 


CONCLUSIONS 

The following are concluded from this study: 

1. A new simplified approximate theoretical method is presented (eq. (20)) for pre- 
dicting the coefficient of thermal contact conductance. It applies to the case of two rough, 
nominally flat, outgassed surfaces in contact and plastically deforming in a high vacuum 
environment. The method does not require the use of a computer program. It does, 
however, require the use of a surface analyzer and a microhardness tester. The analy- 
sis is based on an application of the one- dimensional steady- state Fourier heat conduc- 
tion law to the interface of contact of an idealized ’’button" contact. An effective zone 

of the thermal disturbance is assumed to exist in the vicinity of the interface. The 
zone's minimum length is assumed to be equal to the summation of the average of the 
maximum roughness heights of the two surfaces in contact. These assumptions are 
unique to this method. 

2. The importance of knowing the entire history of events that the test surfaces ex- 
perience during their engagement is emphasized. Thus, meticulous care of the speci- 
men must be taken. Also, a carefully prepared program of experimentation must be 
made and followed to achieve meaningful results. The importance of monotonically in- 
creasing pressure and/or interface temperature has been pointed out. And the necessity 
for uninterrupted testing has been made obvious from the experimental results reported 
here. 

3. An experiment is performed using a tungsten-tungsten specimen tested in a vac- 
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uum environment of 10 to 10 torr. The results of this test are compared with theo- 
retically anticipated results. The agreement is considered good though several deficien- 
cies in the experimental method are pointed out. While the deficiencies in experimental 
method have been pointed out herein, it must also be mentioned that similar deficiencies 
have probably been overlooked in some previous studies. A relation predicting the vari- 
ation of the coefficient of thermal contact conductance during the contact disengagement 
(elastic regime) has been presented. The coefficient may be calculated if the maximum 
pressure at which the contacting surfaces are compressed is known and their mean inter- 
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face temperature remains constant as the pressure is reduced. Also, a relation has 
been presented between the elastic variation of the effective area ratio of contact and 
decreasing pressure. 

4. A new technique for determining the variation of the effective microhardness of 
the interface surface asperities with temperature has been presented. This technique 
has been applied to the tungsten experimental data obtained in the present study. The 
resulting variation of the effective hardness with temperature is shown to agree with 
hardness data in the literature. This technique must be applied after the experimental 
variation of the conductance with pressure has been obtained. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 23, 1969, 

124-09. 
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APPENDIX A 


SYMBOLS 


A 

total cross-sectional area of 

2 

cylindrical model, m 

Pm 

mean pressure on a single asper- 
ity, P m = W'/a;, N/m 2 

A c 

summation of individual areas of 
button contacts (results in total 

Q 

rate of heat transfer, W 


2 

area of actual contact), m 

R 

radius of the cylindrical model, m 

a ; 

2 

contact area of one asperity, 7ra , 
m 2 

r 

radius of curvature of the asperi- 
ties in contact, m 

a 

asperity radius of contact, m 

t 

temperature, K 

C 

constant 

At 

extrapolated temperature differ- 
ence across the interface of 

c 

constant in p m = cY 


contact, K 

d 

diameter of asperity contact, 
d = 2a, m 

At c 

temperature difference across 
the effective zone of the thermal 

d’ 

diameter of planar projection of 
indent resulting from penetra- 


disturbance, K 


tion of a flat surface by a 
spherical indenter 

W 

applied load acting over the total 
nominally flat surface, N 

E 

Young's modulus of elasticity, 
N/m 2 

W' 

load applied to asperity in con- 
tact, N 

h 

coefficient of thermal contact con- 
ductance, kW/(m 2 )(K) 

2 

W o 

upper value of loading in plastic 
range of softer metal (load act- 
ing over total nominally flat 

I 

microhardness, N/m 


surface), N 

k 

coefficient of thermal conductivity, 
W/(m)(K) 

w ; 

upper value of asperity loading in 
plastic range of softer metal 

P 

total load applied to specimen 


(see fig. 2), N 


divided by entire cross - 

2 

sectional area, N/m 

X 

length from interface of button 
contacts along axis of cylinders 

P’ 

load applied to one asperity 


(see fig. 4), m 


divided by the total cross- 
sectional cylindrical area, 
p T = w'/A, N/m 2 

Y 

o 

elastic limit in tension, N/m 
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6 


effective zone of the thermal 
disturbance or average of 
the maximum roughness 
heights of a surface, m 

Subscripts: 

c plane of button contact 


eff effective 

i interface 

o condition of maximum loading 

(see fig. 2) 

1,2, 'l stations or offset in area of 
(T), (2) j contact 
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APPENDIX B 


DISCUSSION OF APPARATUS 

The test apparatus was designed specifically for investigating the thermal contact 
conductance problem. A stainless steel bell jar encloses the test chamber. (Stainless 
steel was used for reasons of safety and cleanliness. ) The jar rests on a base plate and 
on an instrument ring which contains eight vacuum feedthroughs (see fig. 6). An obser- 
vation port located in the bell jar cover allows viewing of the test specimen contact plane 
and some 5.08 centimeters beyond on either side of it. The cover is fastened by bolts 
to the base plate on which an electrically insulated, water cooled, copper heat sink 
rests. A pneumatic loading cylinder (fig. 7) supplies a constant force, denoted by W 
in figure 6, to the cylindrical specimen. The loader is fastened to the top of the bell jar 
by four long stemmed studs which complete a closed force system between the loading 
cylinder and the specimen. 

A small hammer assembly (fig. 7) is used to overcome the static friction of two 
stationary O-ring seals that are located in the top of the bell jar and are in contact with 
the circumference of the loading cylinder heat shield (fig. 6). The hammer is actuated 
intermittently by an electric solenoid through a motor- driven cam and switch. A hand- 

operated hoist located over the facility provides easy removal of the bell jar cover. 
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A vacuum of 10 to 10 torr was maintained in the bell jar during the test. This 
operating condition should be adequate for simulating space conditions. Reference 16 
indicates that (1) the convection of heat conducted to air begins to decrease at approxi- 
mately 70 torr and measurably disappears at 10“* torr, and (2) other investigators have 

-4 

stated that a vacuum in the range of 10 torr is sufficient to negate the conduction ef- 
fect of interstitial fluids. 


Vacuum System 

The vacuum system is composed of a 15. 24-centimeter oil diffusion pump, a me- 
chanical holding pump, high vacuum valves, and the necessary vacuum instrumentation. 
A liquid nitrogen cold trap is located in the high vacuum line between the diffusion pump 
and the bell jar high vacuum valve. The cold trap inhibits the backstreaming of diffusion 
pump oil into the bell jar. The vacuum measuring instrumentation consists of a thermo- 
couple - ionization gage control. 
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Feedthroughs 


Eight feedthrough ports allow flexibility in the system. Two of these are used for 
the passage of water-cooled copper tubes which are connected to the induction coil. An- 
other set of copper cooling lines lead to the copper water-cooled heat sink shown in fig- 
ure 6. A third feedthrough, not shown in figure 6, is used to connect a small liquid- 
nitrogen cryopump to a liquid nitrogen system. The cryopump consists of a 2. 54- 
centimeter-outside- diameter copper finned tube 12. 7 centimeters long which protrudes 
into the bell jar. A fourth feedthrough is used for installation of an ionization gage. All 
of the feedthroughs are sealed to the instrument ring with O- rings. 


Test Specimen and Column 

The following is a consideration of the column of apparatus which includes the test 
specimen. In descending order, this column consists of a pneumatic loader, load cell, 
hammer and solenoid, loading cylinder, induction coil, specimen, heat sink, and base 
plate as shown in figures 6 and 7. 

Loader . - The loader (fig. 7), which is a pneumatic cylinder containing a piston and 
shaft, is used to load the specimen so that data may be obtained as a function of the in- 
terface pressure. Air pressure is applied to the loader by a controller containing a 
solenoid actuated valve which allows air pressure to be applied to either side of the pis- 
ton. The solenoid is actuated by a toggle switch on the control panel (fig. 7). The mag- 

5 

nitude of the pressure is regulated by a pressure regulator with a range 0 to 8. 61x10 

newtons per square meter. The regulator controls the pressure in the piston to within 

about 350 newtons per square meter. This pressure is read out on a direct reading gage 

which has an accuracy of approximately 500 newtons per square meter. 

Load cell . - The load or force on the interface of the specimen is determined by the 

compression of a load cell. The load cell is located between the shaft of the loader and 

the loading cylinder which passes from the outside to the inside of the bell jar. The load 

2 

cell was calibrated by the deadweight method from 0 to 9. 52x10 newtons. The calibra- 

3 

tion was extended to 6. 275x10 newtons through use of a materials testing machine and 
had an accuracy of approximately ±1 newton per square centimeter. The bridge voltage 
of the load cell was maintained at 1.001 volts. This voltage was supplied by an external 
battery supply and rheostat and was measured with the aid of a Rubicon Potentiometer. 
The output voltage of the load cell was read on a nonovoltmeter which had a range from 
0. 1 microvolt to 100 millivolts. Two additional external forces had to be taken into ac- 
count to determine the total true load on the test interface. These forces consisted of 
(1) the weight of the loading cylinder, load cell, and the upper portion of the test speci- 


26 



men (figs. 6 and 7), and (2) the force due to atmospheric pressure acting on the load 

cylinder which has a cross-sectional area of approximately 20. 3 square centimeters. 

o 

The additional weight due to these forces amounted to 2.44x10 newtons. This force re- 
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suited in a pressure of 4. 82x10 newtons per square meter acting on the interface of 
contact in addition to the pressure applied pneumatically. 

Heater coil . - The method used to heat the specimen is of vital concern in a suc- 
cessful test. An inductive heating unit was adapted to the test apparatus and specimen 
arrangement. The inductive coil had 10 turns of hollow copper tubing having an overall 
outside coil diameter of approximately 8.25 centimeters (fig. 6). The coil was water 
cooled from a central manifold which was supplied by city water. The inductive source 
was a 10 080-hertz motor-generator set, which is rated for 220 volts at 136 amperes or 
approximately 30 kilowatts at 0 kilovars. Although the actual input power throughout the 
tests varied from approximately 600 watts to 2. 9 kilowatts, this piece of equipment 
maintained a regulated level of output power within approximately ±3. 5 percent. The 
heat transferred across the test specimen’s interface of contact was varied from 7 to 
400 watts during the test. Since tungsten is nonferrous and has a magnetic permeability 
of approximately 1. 000068 times that of free space, the magnetic field produced by the 
induction coil was unaffected by the test specimen's presence. At the center of the coil's 
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axis, the magnetic field reached a maximum value of approximately 3. 06x10 tesla 
during the test. Though the field induced currents in the vicinity of the specimen's inter- 
face, their effect was negligible. The effect of the induction coil's magnetic field on the 
thermocouples was also negligible. This was because the leads from each thermocouple 
were twisted for some distance away from the zone of the coil and specimen and had no 
loops in them. 

To make use of the existing inductive heating source, the upper portion of the test 
specimen was made long enough so that the contact plane of the interface was 8 centi- 
meters below the lower turn of the inductive heating coil. The specimen's upper portion 
extended approximately 4 centimeters above the coil where it came in contact with a 
1.90-centimeter-diameter sapphire sphere. 

Specimen . - The test specimen was metallurgically prepared from high purity 
pressed and sintered tungsten powder. Its 2. 54-centimeter-diameter cylindrical con- 
figuration was specified because of one- dimensional heat-transfer considerations. The 
specimen's upper portion was 25.4 centimeters long so that it passed completely through 
the 10 turns of the induction heating coil. The lower portion of the specimen was 2. 54 
centimeters long (fig. 6). 

Surface finish . - The lay of the specimen's interface surfaces was produced by 
grinding. The average value of the maximum roughness height across the lay of the 
upper and lower surfaces was 5.35 and 5.51 micrometers, respectively, and their 
arithmetic average roughness heights were 0.538 and 0.943 micrometer, respectively. 
These values were determined after the test was completed. The test was performed 
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with the lay of the top and bottom surfaces approximately parallel. It was found that the 
maximum waviness height (defined in ref. 16) of the surfaces was less than the average 
value of the maximum roughness height. This fact supports the assumption that the 
specimen surfaces were nominally flat. 

Interface temperature and pressure range . - Space environmental conditions of in- 
terest determined the choice of interface temperatures at which the data were obtained. 
Reference 16 indicates that the range of greatest interest extends from 200 to 589 K. 
Because the present investigation only extended upwards from room temperature, the 
values of 394 and 589 K were chosen as minimum and maximum interface temperatures, 
respectively, for the present tungsten-tungsten tests. As noted in figure 8 the actual 
test range extended from approximately 350 to 739 K for the three interface contact 

C £> 

pressures of 4.83x10°, 2.55x10°, and 5.31x10° newtons per square meter. 

Thermometry . - For the design temperatures of 394 to 589 K, a choice of thermo- 
couple materials had to be made. The special type Chromel-Alumel thermocouple wire 
was selected both because of its high millivolt output and its nearly linear voltage varia- 
tion with temperature in the range of interest. These thermocouples have an accuracy 
of ±1.1 K as guaranteed by the manufacturer in the range of temperatures from 273 to 
550 K, and ±3/8 of 1 percent from 550 to 1533 K. The thermocouples were fabricated 
from 30 gage thermocouple wire twisted at its ends and welded. They were then peened 
into holes 0.32 centimeter deep in the surface of the cylinders. The thermocouple 
leads were sheathed in glass insulation. Upon exiting from the holes the leads were 
wound around one-half the specimen's circumference. Several small hanger straps fas- 
tened to the model on one end located the leads on the circumference of the specimen. 

The first straps were located approximately 1. 27 centimeters from the holes to prevent 
interference with the thermal response of the thermocouples in the holes. The output of 
each thermocouple was read on a Rubicon portable potentiometer, which can accurately 
be read to three significant figures and a fourth can be approximated. 

As shown in figure 6, six thermocouples were located on the upper and five on the 
lower portion of the specimen. Also three sets of two thermocouples 180° apart were 
located on the OFHC copper heat meter. A calibration was not made of the thermo- 
couples on the specimen and heat meter. However, figure 13 shows a typical plot of the 
specimen temperature against axial distance, where the readings are generally constant 
with x over moderate temperature ranges. From plots similar to figure 13 the thermo- 
couples were shown to vary within the manufacturers' indicated tolerance for the special 
type of Chromel-Alumel wire used. The data from thermocouple locations 0. 32 centi- 
meter on either side of the interface were not included in the data reduction. This was 
done because they were believed affected by the thermal discontinuity of the interface. 
Reference 28 was used to reduce the raw thermocouple data. 

Heat meter . - As shown in figure 6, a 6. 35 -centimeter oxygen-free high conductivity 
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(OHFC) copper heat meter was located between the lower portion of the test specimen 
and the copper heat sink. The heat-transfer rate through the test specimen was deter- 
mined from the heat meter experimental data. The curve of thermal conductivity against 
temperature used to reduce the copper heat meter data is shown in figure 12. The data 
were obtained from reference 20. 

Heat sink . - As shown in figure 6, the heat meter is screwed into a copper water- 
cooled heat sink. The heat sink contains two internal cooling channels. The maximum 
difference between the inlet and outlet water temperature to the heat sink was estimated 
to be less than 1 K during the test. 

The heat sink rests on a 0. 32 -centimeter thick electric insulator which, in turn, 
rests on the base plate of the bell jar (fig. 6). The electric insulator plus the electrical 
isolation of the water lines leading to the heat sink inhibited an external current system 
from being established between the bell jar and the specimen column. 
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(a) Spherical asperity resting on a 
plane surface. 


^yz^^ZZ777 

(b) Elastic deformation of the contact 
joint. 



(c) Incipient plastic deformation of the 
contact joint. 

Figure 1. 



(d) Plastic deformation of the contact joint. 


- Progressive deformation of a plane surface by spherical asperity. 



Figure 2. - Area of contact as function of load for increasing and decreasing load applica- 
tion to a single asperity indenting a flat plate (ref. 181. 
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^Interface of contact, Tg 



Figure 3. - View of nominally flat contact interface showing the effective zone of the thermal disturbance. 
t 1 >t 14- 


6 = 6 j + 62 (effective zone of 
the thermal disturbance) 



-x -* J »-+x 

x=0 contact plane 


Figure 4. - Sum of areas of contact are shown in form of one button contact and the effective zone of the thermal 
disturbance. 
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— □ — Thermocouple data 

O- Assumed extrapolation of temperature 

data to x = 0 in tx>th cylinders 

Assumed temperature varation across the 

effective zone of the thermal disturbance 


Contact interface 



Figure 5. - Temperature profile along specimen. 
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Figure 7. - High vacuum equipment facility, bell jar, control equipment, and power supply. 







Coefficient of thermal contact conductance, h, kWMrrrMK) 
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Effective area ratio of contact, Ac/A 



Figure 9. - Tungsten -tungsten data of figure 8 replotted as effective area ratio of contact against 
pressure. 
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Effective area ratio of contact, A c /A 



Figure 10. - Selected two and three point tungsten-tungsten data of figure 8 replotted without regard to power 
interruptions. Demonstrates possible misinterpretation of data in plastic regime. 
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Temperature, t, K 

Figure 14. - Thermal conductivity of tungsten. 
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Figure 15. - Microhardness of interface surface of tungsten-tungsten speci- 
men as function of indentor loadings at approximately 294 K. 
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Temperature, t, K 

Figure 16. - Variation of tungsten-tungsten specimen effective microhardness 
as function of temperature computed using data from figures 8 and 15. Other 
tungsten hardness data are from references 25 and 26. 
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